Coulomb gap triptych in a periodic array of metal nanocrystals 
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The Coulomb gap in the single particle density of states (DOS) is a universal consequence of 
electron-electron interaction in disordered systems with localized electron states. Here we show 
that in arrays of monodisperse metallic nanocrystals, there is not one but three identical adjacent 
Coulomb gaps, which together form a structure that we call a "Coulomb gap triptych." We calculate 
the DOS and the conductivity in two- and three-dimensional arrays using a computer simulation. 
Unlike in the conventional Coulomb glass models, in nanocrystal arrays the DOS has a fixed width 
in the limit of large disorder. The Coulomb gap triptych can be studied via tunneling experiments. 



Granular metals and arrays of metallic nanocrystals 
(NCs) represent interesting composite systems, wherein 
the unique properties of individual NCs are combined 
with collective, correlation-driven effects between NCs to 
produce novel material properties [TJ |2] • One of the most 
important of these properties is the electron conductiv- 
ity, which proceeds by electron tunneling, or "hopping" , 
between NCs through the insulating gaps which separate 
them. In relatively dense NC arrays, electron conduction 
can occur both through nearest-neighbor hopping and 
through "cotunneling" of electrons between distant NCs 
via a chain of intermediate virtual states |2] . In the pres- 
ence of some disorder, the latter mechanism dominates at 
low temperatures, where the length of hops grows to op- 
timize the conductivity. This transport mechanism was 
introduced by Mott [3] and is called variable range hop- 
ping (VRH) . When the Coulomb interaction between lo- 
calized electrons is taken into account, it can be shown 
that at sufficiently low temperature VRH conductivity 
obeys the Efros-Shklovskii (ES) law g]: 
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where ctq is a constant (or a weak, power-law function 
of temperature) and Tes is a characteristic temperature. 
Eq. ([T]) has been observed in a number of granular metal 
systems at low temperature (see Refs. [5] and references 
therein). In these systems, as in lightly-doped semicon- 
ductors and other "Coulomb glasses", ES conductivity 
can be seen as the result of a vanishing single-particle 
density of states (DOS) at the Fermi level /i. This van- 
ishing DOS is the consequence of a very general stability 
criterion of the ground state [6] , and it implies that in a 
system of d dimensions the DOS g(E) satisfies 
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Here, Ad is some numerical constant of order unity, E is 
the electron energy relative to the Fermi level, and e is 
the electron charge. Eq. ^ is called the "Coulomb gap." 

In this letter we report an additional striking feature of 
the DOS in periodic arrays of monodisperse metal NCs 
surrounded by random impurity charges. Namely, we 



show that the Coulomb gap at E' = necessarily implies 
the existence of additional, identical Coulomb gaps at 
energies E = ±e^/Co, where Cq is the self-capacitance of 
each NC. This result is shown in Fig. [T] 
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FIG. 1. (Color online) The DOS of a regular array of monodis- 
perse NCs, where E* — E/{e'^ /2Co) is the dimensionless 
single-particle energy and g''{E*) — {e^ D'' /2Co)g{E*) is the 
dimensionless DOS, where D is the NC diameter. Here, the 
results are shown from a computer simulation of a) a 2d square 
lattice and b) a 3d cubic lattice. The shaded area shows 
filled electron states, and the empty area indicates empty 
states. In addition to electron-hole symmetry, the two peaks 
of the DOS have a mirror symmetry across E* = ±1, respec- 
tively (dotted lines). This symmetry creates from the central 
Coulomb gap two additional half-gaps at E* = ±2, resulting 
in a "Coulomb gap triptych." Insets show the DOS near the 
Fermi level i5* = in log-log scale. 



In Fig. [T] one can see that to the right of the Fermi 
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level, at ii^ > 0, the DOS curve g{E) becomes reflected 
across the vertical line E = e^/2Co into the interval 
e^/2Co < E < e^/Co- Similarly, on the left side of the 
Fermi level the portion of the g{E) curve corresponding 
to — e^/2Co < i? < becomes reflected across the ver- 
tical line E ~ — e^/2Co. We refer to this characteristic 
shape as a "Coulomb gap triptych." As we explain below, 
its structure is the result of an additional symmetry in 
the system that arises as a result of the discrete charging 
spectrum of individual NCs. In this way the Coulomb 
gap triptych represents a bridge between the concepts of 
the Coulomb gap and the Coulomb blockade. 

Experimentally, regular arrays of metal NCs can now 
be reliably synthesized with diameter D in the range 3- 
7 nm and with size dispersion less than 5% [TJ |2l [S]. 
For such small NCs, the self-capacitance Co is also small: 
Co = kD/2, where k is the effective dielectric constant of 
the array, given approximately by the Maxwell-Garnett 
formula fT, '5^ . Correspondingly, the Coulomb self-energy 
q^/2Co of an NC with charge q plays a large and impor- 
tant role in electron transport. To see this, one can imag- 
ine a hypothetical NC array with no disorder. In such an 
array, in the ground state all NCs are neutral and electron 
conduction requires the thermal excitation of positive- 
negative NC pairs. Thus, the conductivity is activated 
with an activation energy e^/2Co. For nanometer-sized 
NCs, this activation energy can easily exceed the thermal 
energy ksT. 

In the presence of some finite charge disorder, however, 
the fluctuating Coulomb potential can produce charging 
of NCs in the ground state and thus lead to a Coulomb 
gap in the DOS and to ES conductivity. To show how 
this happens, in this letter we adopt the following simpli- 
fied model. We assume that identical, spherical, metallic 
NCs reside in a regular d-dimensional square lattice with 
lattice constant D' , and that impurity charges ±e are 
embedded in the insulator (oxide) between NCs. Such 
impurity charges can be thought to effectively create a 
fractional donor charge Qi that resides on each NC i, for 
reasons that are explained below. The net charge of the 
NC can then be written as qi = Qi — erii, where rii is the 
integer number of electrons that reside on the NC rela- 
tive to its neutral state {rii can be positive or negative). 
Given this model, the Hamiltonian for the system is 

g - E C-,\Q.-en^m-^n,). (3) 
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Here, the first term describes the Coulomb self-energy of 
each NC and the second term describes the interaction 
between charged NCs. The coefficient C^j^ is the inverse 
of the matrix of electrostatic induction Cij . This Hamil- 
tonian has been also been proposed as a model for arrays 
of large semiconductor NCs [9] . 

Because of the presence of the impurity charges, elec- 
trons become redistributed among NCs from their neu- 



tral state in order to screen the disorder Coulomb po- 
tential. In order to calculate the DOS and conductivity 
we first attempt to find numerically the set of electron 
occupation numbers {rii} that minimizes the Hamilto- 
nian. In the numerical simulations that we describe be- 
low, we make the approximations that Cq — kD/2 and 
C^j^ = l/KVij. These approximations do not effect our 
main conclusions, as we explain below. 

The model of fractional donor charges Qi was first put 
forward in Ref. [TU]; here its justification is briefly re- 
peated. When an impurity charge, say with charge +e, 
is located close to the point of contact between two NCs, 
labeled A and B, it induces negative image charges —qA 
and —qB in the surfaces of NCs A and B, respectively. 
This is shown schematically in Fig. [2] In order to main- 
tain overall neutrality of the NCs, an equal and opposite 
image charge appears at the center of each NC: +qA and 
+qB- (These "image charges at the center" represent a 
uniform electronic charge at the NC surface.) The values 
of qA and qB are such that together the image charges 
—qA and —qB neutralize the donor charge: qA + qs = e. 
Their respective magnitudes are determined by the dis- 
tance between the impurity and each NC surface. For ex- 
ample, if the impurity sits exactly along the line connect- 
ing the centers of NCs A and B and if the gap w = D' — D 
between NCs satisfies w <^ D, then qAXB = qnXA, where 
XA and xb are the distances between the impurity and 
the surface of NCs A and B, respectively. Since the im- 
purity charges and the image charges —qA, —qs together 
form a compact, neutral arrangement, the net effect of 
the impurity charge is to produce "fractionalized" donor 
charges, such that +qA is relayed to the center of NC A 
and +qB is relayed to the center of B. 




FIG. 2. (Color online) A schematic depiction of the fraction- 
alization of a charged impurity (small black circle) between 
NCs (large gray circles). The positive impurity induces neg- 
ative image charges (white circles) in nearby metal surfaces 
and is effectively neutralized, while equal and opposite posi- 
tive images are conveyed to the center of the NC (x's). 

In this way, each NC i can be said to have a frac- 
tional donor charge Qi, which is equal to the sum total 
of the fractionalized charges donated by individual impu- 
rities around it. In the limit where there are very many 
impurity charges surrounding each NC, one can think 
that the random variable Qi is Gaussian-distributed with 
some standard deviation larger than e. In fact, how- 
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ever, in such cases one can effectively adopt a much sim- 
pler model, in which the value of Qi is chosen randomly 
from the uniform distribution Qi e [— e/2,+e/2]. To see 
why this model is valid, consider that each NC minimizes 
its Coulomb self-energy by minimizing the magnitude of 
its net charge, \Qi — eni\. Since can take any inte- 
ger value, it is generally true that in the ground state 
— e/2 < Qi — erii < e/2. In other words, each NC can 
effectively adjust to the presence of an arbitrarily strong 
charge disorder by changing its electron number (say, 
by drawing electrons from the voltage source) so that its 
net charge acquires a magnitude smaller than e/2. This 
has important implications for the disorder-dependence 
of conductivity, as we show below. 

Given the ground state configuration for a particular 
system, defined by the set of electron occupation numbers 
{rii}, one can determine the energy of the highest filled 
electron level, and the lowest empty electron level, 

E^''\ at each NC i. Specifically: 
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These energies are defined so that the Fermi level /i = 0, 
and in the ground state e\^^ < and e'^''^ > for all 
i. The single particle DOS g{E) is defined by making 
a histogram of the energy values E^^'^^ and e'^'^^ . Higher 
and lower electron energy states are ignored in this work, 
as they play no role in conductivity at ksT <^ /Cq. 

In order to evaluate numerically the DOS, we use a 
computer simulation to search for the ground state ar- 
rangement of electrons, {n^}, in a finite array of NCs. 
For simplicity, we set the lattice constant D' — D; this 
corresponds to the limit where the gap w between NCs 
is very thin while the tunneling transparency of the bar- 
rier between them remains much smaller than unity. In 
our simulation we search for the ground state by looping 
over all NC pairs i,j and attempting to move one elec- 
tron from i to j. If the move lowers the Hamiltonian H, 
then it is accepted, otherwise it is rejected. Equivalently, 
one can say that for all i, j we check that the ES ground 
state criterion is satisfied: 
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eP - e^CZ^ > 0. 
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It should be noted that this procedure does not in general 
find the exact ground state, but only a "pseudo-ground 
state" that is stable with respect to single-electron trans- 
fers. In principle, the system energy can be lowered fur- 
ther by some simultaneous multi-electron transfers. Such 
processes are generally seen to have only a relatively 
weak effect on the DOS [11] [T2] that slightly deepens 
the Coulomb gap near the Fermi level. 



The resulting DOS is shown in Fig. [T^ for a two- 
dimensional (2d) simulated system of size 100 x 100 lat- 
tice sites and in Fig. [TJd for a three-dimensional (3d) sys- 
tem of size 25 X 25 X 25. Electron energies are plot- 
ted in the dimensionless form E* = E / {e^ /2Cq) and 
the DOS is plotted in the dimensionless form g*[E*) = 
{e^D'^ /2Co)g{E*). The insets to these figures show a 
log-log plot of the DOS near E = Q, which suggests that 
in 2d the DOS follows g2d{E) oc E^-^ at small energies 
and in 3d gsdiE) cx E'^ '^. These exponents are some- 
what larger than the theoretical ones given in Eq. ([T]), so 
that apparently the ES bound is not saturated. This 
is similar to what happens in the Efros model of the 
Coulomb glass [13] at disorder strength A = 1 ^Tj. The 
results of Fig. [T] are generated using a uniform distribu- 
tion Qi g [—e/2, e/2] for the fractional charge. If one 
instead takes Qi to be Gaussian-distributed with a stan- 
dard deviation < 3e, the resulting DOS is everywhere 
equal to that of Fig. [T] to within 0.6%. 

Fig. [T] also highlights the striking additional symme- 
try in the DOS in both 2d and 3d, as compared to the 
DOS in the conventional Coulomb glass problem [51 lllj. 
Namely, each peak in the DOS is symmetric with respect 
to reflections about E* = ±1, so that the DOS has iden- 
tical, repeated Coulomb gaps at E* — ±2. The origin 
of these additional Coulomb gaps can be understood by 
noting a particular symmetry in the Hamiltonian that is 
reflected in the filled and empty state energies. El" and 
E^'^K Namely, by subtracting Eqs. Q and ^ one can 
show that 



E*'^'^ = e: 
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for all i. Thus, all NCs contribute to the DOS two energy 
levels - one filled, one empty - separated by e^/Cg. This 
implies that as the density of states collapses at E very 
close to zero (the Coulomb gap), the density of states 
must also collapse as E* approaches ±2 in identical fash- 
ion. That is, the ES stability criterion of Eq. ([6| places 
constraints both on the DOS near E = and on the DOS 
near E = ±e^/Co. 

One can also note that states with E*^-^^ < —2 or 
-E*^^'' > 2 are prohibited, since by Eq. hi these would 



imply that some NC has E^'"'^ < or eJ^^ > 0. Thus, 
g{E) is strictly zero at \E*\ > 2. This is a markedly 
different situation than in the conventional Efros model 
|13j . where the width of the DOS reflects the character- 
istic strength of the disorder. In the present problem, 
for large enough disorder the DOS has a saturated width 
e^/Cp. This saturation occurs because the number of 
electrons n at each site can adjust to screen an arbitrar- 
ily large Coulomb disorder. Thus, one can expect that at 
large disorder the conductivity also becomes independent 
of disorder strength. 

In order to evaluate the conductivity directly, we em- 
ploy the approach of the Miller- Abrahams network [14j . 
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in which each pair ij of NCs is said to be connected 
by some equivalent resistance Rij . The value of Rij in- 
creases exponentially with the distance rij between NCs 
and the activation energy AEij required for electron hop- 
ping between i and j according to Rij (x exp[2rij/^ + 
AEij /ksT]^ where ^ is the electron localization length 
[TU] and the value of AE'ij is determined by the ground 
state energies {E^^^} and {e\'^^} [5]. The resistance of 
the system as a whole can be found using a percolation 
approach. Specifically, we find the minimum value Rc 
such that if all resistances Rij with Rij < Rc are left in- 
tact while others are eliminated (replaced with i? = oo), 
then there exists a percolation pathway connecting op- 
posite faces of the simulation volume. The conductivity 
of the system a is equated with 1 / {R^D'^-'^) . 

Our results for the conductivity are shown in Fig. [Sj 
plotted as a function of the dimensionless temperature 
T* = 4DCokBT/{e^^) raised to the power -1/2. The re- 
sults indicate that the conductivity is well-described by 
the ES law of Eq. ([T]) at relatively small temperatures 
T* < 1, both in 2d and 3d ^5J. This behavior is consis- 
tent with the prominent Coulomb gaps seen in Fig. [l] In 
both 2d and 3d, replacing the uniform distribution of Qi 
with a distribution with larger variance — for example, 
by taking Qi as the sum of three or more independent 
fractional charges — did not affect the conductivity to 
within our numerical accuracy. This insensitivity to the 
disorder strength stands in contrast to the Efros model 
[13], where large disorder widens the DOS, so that ES 
conductivity exists only when the temperature is suffi- 
ciently small that electron hops are confined to within the 
parametrically narrow window of energies in which g{E) 
is constrained by the Coulomb gap On the contrary, 
in arrays of monodisperse metallic NCs the DOS becomes 
essentially independent of disorder strength, so that even 
at large disorder the Coulomb gap plays a prominent role 
and the conductivity follows the ES law. 

The triptych structure of the DOS should have observ- 
able consequences for a number of experiments on metal 
NC arrays. It is possible, for example, that the DOS 
can be probed directly by tunneling experiments, simi- 
lar to the ones that have directly observed the Coulomb 
gap in doped semiconductors [TB]. For systems with a 
finite dispersion 5C in the NC self-capacitance, the re- 
peated Coulomb gaps will be smeared over some finite 
energy interval rather than collapsing to zero exactly 
at E* = ±2. One can simulate this behavior numeri- 
cally by adding a stochastic spatial variation to Cq. Our 
simulations suggest that for root mean square deviation 
SC < Co, g{E* = ±2)/g{E* = ±1) w 3{SC/Cof. This 
implies that for a system with 5% dispersion in the NC 
diameter, the collapse of the DOS at E* = ±2 is com- 
plete to within 1%, and the resulting g*{E*) curve would 
not be distinguishable from that of Fig. [T] if added to the 
plot. 
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FIG. 3. (Color online) The temperature dependence of the 
conductivity in a) 2d and b) 3d. In both cases, the con- 
ductivity follows the ES law [Eq. ([T])] at small temperatures, 
T* ^ 1, as shown by the dashed lines. 
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